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Abstract 

The reasons for studying cosmic-ray transport theory are 
summarized and the fundamentally three-dimensional nature of the 
process is pointed out. Observations in the ecliptic plane cannot 
unambiguously test transport theories since the solutions to the 
transport equations depend critically on boundary conditions and 
variation of parameters such as diffusion tensor out of the ecliptic. 
Sample calculations are shown which illustrate the problem. It is 
concluded that out-of-the-eci iptic observations are essential to 
further test transport theory. 
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I. Introduction 

The study of cosmic-ray transport theory has an old and vener- 
able history, for It is clear that In order to understand anything 
in cosmic-ray astrophysics we must first understand transport. The 
observed near-isotropy of cosmic rays was early recognized to Imply 
that the orbits of the particles have been severely distorted in 
their motion through space. Since the work of Fermi (19**9) It has 
been understood that the motion must be treated statistically, since 
the plasmas and magnetic fields through which the cosmic rays move 
are Irregular and turbulent. This basic fact leads to the view that 
the spatial motion of cosmic rays is to a very good first approxima- 
tion a random walk in three dimensions. 

The solar wind provides an excellent local laboratory for testing 
our ideas of cosmic-ray transport by comparing in situ observations 
with theory. A comprehensive theory has been developed which has 
had reasonable success in explaining the various cosmic-ray phenomena 
in the solar wind (see reviews by Jokipii, 1971 and Volk, 1975). 

Before going on to discuss this theory in detail, it is useful to 
emphasize that irrespective of the detailed theory, cosmic-ray trans- 
port theory is fundamentally three-dimensional and that this problem 
is much more severe than for some other aspects of the physics of 
the interplanetary medium. The general picture is illustrated in 
figure I for particles originating in the galaxy or at the sun. It 
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one sees cosmic rays which have sampled only the ecliptic piane. 

This conclusion depends only on the turbulent nature of the inter- 
planetary plasma and the consequent stochastic motion of the cosmic- 
ray particles. 

One can, of course, dream up experiments which could be used to 
study cosmic-ray transport locally. For example, a useful experi- 
ment would be to inject a smali beam of "tagged" cosmic rays at one 
point and measure them a distance -0.1 a.u. away. But such experi- 
ments are clearly not possible in the near future. We must be 
content with v.'nat nature provides, and work with the full three- 
dimensional problem. 

The preceding discussion hjs intentionally been as general as 
possible in order to emphasize the fundamental nature 
O' the conclusions. In the next section, the current detailed theory 
of cosmic-ray transport is discussed, then some quantitative calcu- 
lations are presented which illustrate the effects of uncertainties 
in parameter outside of the ecliptic. 
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II. Current Transport Theory 


The current theory of cosmic-ray transport has been comprehen- 
sively discussed In an earlier review by the author (Joklplt, 1970 • 
The theory works at two levels, which could be termed the "macro- 
scopic" and the "microscopic". This Is Illustrated schematically 
In figure 2 which shows the particle being scattered randomly by the 
magnetic fluctuations. Hence the particle pitch angle 9 undergoes a 
random walk. The details of this scattering process and Its relation 
to the detailed microstructure of the magnetic fluctuations are 
studied In the microscopic theory, whereas the resulting spatial 
diffusion is considered in macroscopic theory. Thus, for example, 
in the simplest form of quasi linear theory, where only planar magnetic 
fluctuations with wave vector parallel to the average field are con- 
sidered, the pitch-angle scattering Is characterized by the Fokker- 
Planck coefficient 


<A0^> 

At 




( 1 ) 


v/here ■ qB^/mc Is the particle cyclotron frequency in the average 
field Bjj, y = COS0, w » particle speed, r^ ■ w/Wj,, and Pj_(k) is the 
spatial power spectrum of the magnetic fluctuations as a function of 
wavenumber k (see, e.g. , Jokipii 1966, 1971). Other more complex 
expressions result if other magnetic fluctuation configurations are 
assumed. Similar expressions result for the other Fokker-Planck 
coefficients such as <Ax^>/At, etc. Hence, if x is a direction 
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normal to we have 


<Ax^> 

At 


■^P, (k«0) . 
B 2 


( 2 ) 


The above expressions are examples of microscopic transport theory. 

However, in most situations the cosmic-ray angular distribution 

is driven essentially isotropic by the scattering, so that the 

transport must be described by the diffusion approximation. If 

U(r,t,T) Is the particle density averaged over pitch angle as a 

function of position r, time t rnd energy T, then the diffusion 

equation reads in the rest frame (Parker 1965, Gleeson and Axford 1967, 
Jokipii and Parker 1970) 
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with an associated flux of particles in the rest frame 
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where <jj is the cosmic-ray diffusion tensor, is the solar wind 
velocity and o(T) » (2mc^ + T)/(mc^ + T). The associated anisotropy 
1^1 “ 3|Fj/wU. In a coordinate system with the z-direction oriented 
along the average magnetic field B^, k.j may be written 
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where the parallel dtffuston coefficient K|, , the perpendicular 
diffusion coefficient Kj_and the antisymmetric diffusion coefficient 
can be written 





( 6 ) 

(7) 


The diffusion equation (3) embodies the macroscopic theory which is 
connected to the microscopic theory by equations (5) " (7). At 
present, comparison of transport theory with observation must be 
done through the use of macroscopic theory and then working back to 
the microscopic theory. We are unable to measure <A^/At, etc. 
d I rectly. 

With regard to the problem of out-of-the-ecl iptic exploration, 
it Is clear that the solution to equation (3) depends on a knowledge 
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of K|j and througnout the modulating region and on the value of 
U on the boundary. There appear to be no purely locally measurable 
properties of the solution which can be used in checking transport 
theory. 
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III. MlustratWe Calculations 


To Mlustrata the rather large uncertainties Introduced by lack 
of knowledge of parameters out of the ecliptic, I present In this 
section a summary of some analytical calculations published elsewhere 
(Owens and Jokipll, 1971}. This problem has also been considered by 
Sarabhal and Subramanlan ( 1966 ) and LlettI and Quenby (I 968 ). Con- 
sider th'i 11-^year solar cycle modulation of galactic cosmic rays by 
a solar wind which Is not spherically symmetric. The 11-year solar 
cycle variations are usual?'’ regarded as slow enough that the time 
derivative In equation ( 3 ) can be neglected in solving for U. In 
this calculation the various parameters In the transport equation 
are allowed to vary with heliographic latitude and the effects of 
the variation on the density U and the flux F in the solar equatorial 
plane are considered. 

The following forms for the various parameters were assumed: 

The various components of the diffusion tensor are Independent of 
energy T and are proportional to heliocentric radius r out to some 
boundary r ■ 0 where -► «. The cosmic-ray density U^r,T) is 
assumed t:- take on a given Interstellar value U (T) - AT“^*5 at r « 0. 

09 

Finally, we follow the us;;^1 practice of circumventing the 
problem of proper boundary conditions at the Sun by requiring the 
solution for U to be finite at the origin. This will be shown to 
lead to negligible error in the present case. See Jokipii (1971) 
for a more complete discussion of this problem. 
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Equation (3a) for U(r,0,$,T) becomes, tn spherical polar 
coorJi nates 
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( 8 ) 


in what follows, we consider only mocels which are axisymnetric 
(independent of ^). If 0 * 0 is the axis of solar rotation, this 
corresponds to a lati tude-dependent solar wind. Effects of magnetic 
sectors could be represented by taking 0=0 along a sector. The 
two terms containing cancel from the density equation (8) if 
3<^/33 = 3x^/3$ = 0, as in our axi symmetric model. Thus the terms 
in :<^ do not contribute, and equation (8) becomes 
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Since x„and Xj^ are taken to be independent of energy, equation (9) 
separates, and one writes 


U(r,Q,T) = S(r,0)5(T). 


(10) 
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Upon neglecting the small energy dependence of a(T), one finds 
immediately that S must satisfy 4- 



with 

0(T) - T(3q/2a)-l . (12) 

The separation constant q must be chosen to agree with 
U^(T) « T“2.5^ }p which case we choose q » -1. 

A number of different forms for the variation of the parameters 
V* ''•* chosen. Define 

“ V^[l+5(0)J, Kj|(0,r) » K^[l+e(0)]r 

and “ Hfii • ( 13 ) 

The problem may be solved in terms of a series of Legendre 
polynomials and the general solution is given by Oweis and Jokipii 
( 1971 ) •. An illustrative velocity variation is given by 
V « V^[l + .BOP^CcosG)] where is the average velocity and P 2 is 

the Legendre polynomial of order 2. This velocity variation is 
illustrated In figure 3- Some typical results of the calculations 
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are illustrated iu figuras 4, 5> and 6. - it is clear that many of 
the parameters observed in the solar equatorial plane can be sub- 
stantially changed by varying the solar-wind parameters out of the 
solar ecliptic plane. Of particular note is the fact that the radial 
anisotropy observed in the equatorial plane is extremely sensitive 
to parameters outside the equatorial plane. As shown in figure 6, 
even the sign of the anisotropy can be changed by relatively small 
velocity variations. 
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IV. Conclusions 

may conclude from the above discussion that fundamental 
ambiguities in testing cosmic-ray transport can be removed by 
carrying out measur^ients out of the ecliptic. The out-of-the-eci iptic 
measurements most necessary to study cosmic-ray propagation are: 


a) Measurements of the flux of cosmic rays as a function 
of solid angle and energy with as much resolution as 
possible. The anisotropy measurements are easier to 
carry out on a spinning spacecraft. 

b) Simultaneous measurements on board the same spacecraft 
of the plasma and magnetic field. These measurements 
should be spaced in time so that good time-series 
analyses (power spectra, etc.) can be obtained. 

c) The spacecraft should go as far out of the ec 1 iptic 
as possible to insure that any variations with helio- 
graphic latitude will be seen. 

d) Good base-line measurements at Earth or in the eclip- 
tic plane should be obtained simultaneously. 


It is not crucial that these measurements be carried out at 
constant heliocentric radius, although this might aid interpretation. 
It appears that measuring protons and possibly electrons will be 
adequate and it is better to optimize measurements for one species 
rather than compromise these in order to study composition. 
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Figure Captions 


Fig. 1. 


Fig. 2. 

Fig. 3. 
Fig. h. 


Fig. 5. 


Schematic illustration in a meridian plane of typical 
trajectories of solar cosmic-rays (dashed line) and 
galactic cosmic rays (solid line) which reach a detector 
at I a.u. The fundamentally 3'dimensional character of 
the motion is apparent. 

illustration of the "scattering" of a cosmic-ray particle 

by a magnetic irregularity. If the region over which A0 

is correlated is small, a Fokker-Pianck equation results. 

illustration of the velocity profile used in some of the 

calculations. 3 is the angle relative to the z-axis. 

Comparison of linearized and exact solutions for the case 

V = V^£l + 0.30 p 2 (cos 0)], Kj. = Kj^r, p = 3. The density 

U(r,9) is given for both solutions as a function of angle 0 

r = 0.2 D in terms of = U(r = 0,0). Straight line shows 

the corresponding result for an isotropic wind V = V^. 

The linearized calculation keeps terms only to first order 

in 6 as defined in equation (l») . From Owens and Jokipii (1971). 

Total cosmic-ray flux as a function of r,0, and with Fx 

9 

suppressed. The pattern is azimuthal ly symmetric and even 

about the equator, (a) The parameters V(o) = Vp[l + 0.30 P 2 (cos 0)], 

Kj^(r,o) = x^r and p = j v/ere ubed. The position of Earth 

for 0=5 a.u. is indicated. Of particular interest is the 
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virtual source of particles in the equatorial plane at 
r “ 0.5 D. (b) Same as (a) except that V(0) = Vjj[l + 0.46 P^Ccos 0)]. 
Note that the flux Is much larger and virtual source has 
moved out to r ~ 0.9 0. From Owens and Jokipll (1971). 

Fig. 6. Anisotropy In the equatorial plane at r ■ 0.2 0 for 
V(0) a + V 2 P 2 (cos O) , as a function of V^/V^. 

Note that In our model the anisotropy Is radial. From 
Owens and Jokipll (1971) • 
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